Quasinonlocal coupling of nonlocal diffusions

We developed a new self-adjoint, consistent, and stable coupling strategy for nonlocal diffusion models, inspired by the quasinonlocal atomistic-to-continuum method for crystalline solids. The proposed coupling model is coercive with respect to the energy norms induced by the nonlocal diffusion kernels as well as the $L^2$ norm, and it satisfies the maximum principle. A finite difference approximation is used to discretize the coupled system, which inherits the property from the continuous formulation. Furthermore, we design a numerical example which shows the discrepancy between the fully nonlocal and fully local diffusions, whereas the result of the coupled diffusion agrees with that of the fully nonlocal diffusion.Comment: 28 pages, 3 figures, ams.or

A quasinonlocal coupling method for nonlocal and local diffusion models

In this paper, we extend the idea of "geometric reconstruction" to couple a nonlocal diffusion model directly with the classical local diffusion in one dimensional space. This new coupling framework removes interfacial inconsistency, ensures the flux balance, and satisfies energy conservation as well as the maximum principle, whereas none of existing coupling methods for nonlocal-to-local coupling satisfies all of these properties. We establish the well-posedness and provide the stability analysis of the coupling method. We investigate the difference to the local limiting problem in terms of the nonlocal interaction range. Furthermore, we propose a first order finite difference numerical discretization and perform several numerical tests to confirm the theoretical findings. In particular, we show that the resulting numerical result is free of artifacts near the boundary of the domain where a classical local boundary condition is used, together with a coupled fully nonlocal model in the interior of the domain

NySALT: Nystr\"{o}m-type inference-based schemes adaptive to large time-stepping

Large time-stepping is important for efficient long-time simulations of deterministic and stochastic Hamiltonian dynamical systems. Conventional structure-preserving integrators, while being successful for generic systems, have limited tolerance to time step size due to stability and accuracy constraints. We propose to use data to innovate classical integrators so that they can be adaptive to large time-stepping and are tailored to each specific system. In particular, we introduce NySALT, Nystr\"{o}m-type inference-based schemes adaptive to large time-stepping. The NySALT has optimal parameters for each time step learnt from data by minimizing the one-step prediction error. Thus, it is tailored for each time step size and the specific system to achieve optimal performance and tolerate large time-stepping in an adaptive fashion. We prove and numerically verify the convergence of the estimators as data size increases. Furthermore, analysis and numerical tests on the deterministic and stochastic Fermi-Pasta-Ulam (FPU) models show that NySALT enlarges the maximal admissible step size of linear stability, and quadruples the time step size of the St\"{o}rmer--Verlet and the BAOAB when maintaining similar levels of accuracy.Comment: 26 pages, 7 figure

Multi-decadal trends in global terrestrial evapotranspiration and its components

Evapotranspiration (ET) is the process by which liquid water becomes water vapor and energetically this accounts for much of incoming solar radiation. If this ET did not occur temperatures would be higher, so understanding ET trends is crucial to predict future temperatures. Recent studies have reported prolonged declines in ET in recent decades, although these declines may relate to climate variability. Here, we used a well-validated diagnostic model to estimate daily ET during 1981–2012, and its three components: transpiration from vegetation (Et), direct evaporation from the soil (Es) and vaporization of intercepted rainfall from vegetation (Ei). During this period, ET over land has increased significantly (p < 0.01), caused by increases in Et and Ei, which are partially counteracted by Es decreasing. These contrasting trends are primarily driven by increases in vegetation leaf area index, dominated by greening. The overall increase in Et over land is about twofold of the decrease in Es. These opposing trends are not simulated by most Coupled Model Intercomparison Project phase 5 (CMIP5) models, and highlight the importance of realistically representing vegetation changes in earth system models for predicting future changes in the energy and water cycle

Monolayer Excitonic Laser

Recently, two-dimensional (2D) materials have opened a new paradigm for fundamental physics explorations and device applications. Unlike gapless graphene, monolayer transition metal dichalcogenide (TMDC) has new optical functionalities for next generation ultra-compact electronic and opto-electronic devices. When TMDC crystals are thinned down to monolayers, they undergo an indirect to direct bandgap transition, making it an outstanding 2D semiconductor. Unique electron valley degree of freedom, strong light matter interactions and excitonic effects were observed. Enhancement of spontaneous emission has been reported on TMDC monolayers integrated with photonic crystal and distributed Bragg reflector microcavities. However, the coherent light emission from 2D monolayer TMDC has not been demonstrated, mainly due to that an atomic membrane has limited material gain volume and is lack of optical mode confinement. Here, we report the first realization of 2D excitonic laser by embedding monolayer tungsten disulfide (WS2) in a microdisk resonator. Using a whispering gallery mode (WGM) resonator with a high quality factor and optical confinement, we observed bright excitonic lasing in visible wavelength. The Si3N4/WS2/HSQ sandwich configuration provides a strong feedback and mode overlap with monolayer gain. This demonstration of 2D excitonic laser marks a major step towards 2D on-chip optoelectronics for high performance optical communication and computing applications.Comment: 15 pages, 4 figure

Evaluating Alternative Ebullition Models for Predicting Peatland Methane Emission and Its Pathways via Data–Model Fusion

Understanding the dynamics of peatland methane (CH4) emissions and quantifying sources of uncertainty in estimating peatland CH4 emissions are critical for mitigating climate change. The relative contributions of CH4 emission pathways through ebullition, plant-mediated transport, and diffusion, together with their different transport rates and vulnerability to oxidation, determine the quantity of CH4 to be oxidized before leaving the soil. Notwithstanding their importance, the relative contributions of the emission pathways are highly uncertain. In particular, the ebullition process is more uncertain and can lead to large uncertainties in modeled CH4 emissions. To improve model simulations of CH4 emission and its pathways, we evaluated two model structures: (1) the ebullition bubble growth volume threshold approach (EBG) and (2) the modified ebullition concentration threshold approach (ECT) using CH4 flux and concentration data collected in a peatland in northern Minnesota, USA. When model parameters were constrained using observed CH4 fluxes, the CH4 emissions simulated by the EBG approach (RMSE = 0.53) had a better agreement with observations than the ECT approach (RMSE = 0.61). Further, the EBG approach simulated a smaller contribution from ebullition but more frequent ebullition events than the ECT approach. The EBG approach yielded greatly improved simulations of pore water CH4 concentrations, especially in the deep soil layers, compared to the ECT approach. When constraining the EBG model with both CH4 flux and concentration data in model–data fusion, uncertainty of the modeled CH4 concentration profiles was reduced by 78 % to 86 % in comparison to constraints based on CH4 flux data alone. The improved model capability was attributed to the well-constrained parameters regulating the CH4 production and emission pathways. Our results suggest that the EBG modeling approach better characterizes CH4 emission and underlying mechanisms. Moreover, to achieve the best model results both CH4 flux and concentration data are required to constrain model parameterization

Disentangling land model uncertainty via Matrix-based Ensemble Model Inter-comparison Platform (MEMIP)

Background Large uncertainty in modeling land carbon (C) uptake heavily impedes the accurate prediction of the global C budget. Identifying the uncertainty sources among models is crucial for model improvement yet has been difficult due to multiple feedbacks within Earth System Models (ESMs). Here we present a Matrix-based Ensemble Model Inter-comparison Platform (MEMIP) under a unified model traceability framework to evaluate multiple soil organic carbon (SOC) models. Using the MEMIP, we analyzed how the vertically resolved soil biogeochemistry structure influences SOC prediction in two soil organic matter (SOM) models. By comparing the model outputs from the C-only and CN modes, the SOC differences contributed by individual processes and N feedback between vegetation and soil were explicitly disentangled. Results Results showed that the multi-layer models with a vertically resolved structure predicted significantly higher SOC than the single layer models over the historical simulation (1900–2000). The SOC difference between the multi-layer models was remarkably higher than between the single-layer models. Traceability analysis indicated that over 80% of the SOC increase in the multi-layer models was contributed by the incorporation of depth-related processes, while SOC differences were similarly contributed by the processes and N feedback between models with the same soil depth representation. Conclusions The output suggested that feedback is a non-negligible contributor to the inter-model difference of SOC prediction, especially between models with similar process representation. Further analysis with TRENDY v7 and more extensive MEMIP outputs illustrated the potential important role of multi-layer structure to enlarge the current ensemble spread and the necessity of more detail model decomposition to fully disentangle inter-model differences. We stressed the importance of analyzing ensemble outputs from the fundamental model structures, and holding a holistic view in understanding the ensemble uncertainty

Evaluating Alternative Ebullition Models for Predicting Peatland Methane Emission and Its Pathways via Data–Model Fusion

Understanding the dynamics of peatland methane (CH4) emissions and quantifying sources of uncertainty in estimating peatland CH4 emissions are critical for mitigating climate change. The relative contributions of CH4 emission pathways through ebullition, plant-mediated transport, and diffusion, together with their different transport rates and vulnerability to oxidation, determine the quantity of CH4 to be oxidized before leaving the soil. Notwithstanding their importance, the relative contributions of the emission pathways are highly uncertain. In particular, the ebullition process is more uncertain and can lead to large uncertainties in modeled CH4 emissions. To improve model simulations of CH4 emission and its pathways, we evaluated two model structures: (1) the ebullition bubble growth volume threshold approach (EBG) and (2) the modified ebullition concentration threshold approach (ECT) using CH4 flux and concentration data collected in a peatland in northern Minnesota, USA. When model parameters were constrained using observed CH4 fluxes, the CH4 emissions simulated by the EBG approach (RMSE = 0.53) had a better agreement with observations than the ECT approach (RMSE = 0.61). Further, the EBG approach simulated a smaller contribution from ebullition but more frequent ebullition events than the ECT approach. The EBG approach yielded greatly improved simulations of pore water CH4 concentrations, especially in the deep soil layers, compared to the ECT approach. When constraining the EBG model with both CH4 flux and concentration data in model–data fusion, uncertainty of the modeled CH4 concentration profiles was reduced by 78 % to 86 % in comparison to constraints based on CH4 flux data alone. The improved model capability was attributed to the well-constrained parameters regulating the CH4 production and emission pathways. Our results suggest that the EBG modeling approach better characterizes CH4 emission and underlying mechanisms. Moreover, to achieve the best model results both CH4 flux and concentration data are required to constrain model parameterization

Evaluating alternative ebullition models for predicting peatland methane emission and its pathways via data–model fusion

Understanding the dynamics of peatland methane (CH4) emissions and quantifying sources of uncertainty in estimating peatland CH4 emissions are critical for mitigating climate change. The relative contributions of CH4 emission pathways through ebullition, plant-mediated transport, and diffusion, together with their different transport rates and vulnerability to oxidation, determine the quantity of CH4 to be oxidized before leaving the soil. Notwithstanding their importance, the relative contributions of the emission pathways are highly uncertain. In particular, the ebullition process is more uncertain and can lead to large uncertainties in modeled CH4 emissions. To improve model simulations of CH4 emission and its pathways, we evaluated two model structures: (1) the ebullition bubble growth volume threshold approach (EBG) and (2) the modified ebullition concentration threshold approach (ECT) using CH4 flux and concentration data collected in a peatland in northern Minnesota, USA. When model parameters were constrained using observed CH4 fluxes, the CH4 emissions simulated by the EBG approach (RMSE = 0.53) had a better agreement with observations than the ECT approach (RMSE = 0.61). Further, the EBG approach simulated a smaller contribution from ebullition but more frequent ebullition events than the ECT approach. The EBG approach yielded greatly improved simulations of pore water CH4 concentrations, especially in the deep soil layers, compared to the ECT approach. When constraining the EBG model with both CH4 flux and concentration data in model–data fusion, uncertainty of the modeled CH4 concentration profiles was reduced by 78 % to 86 % in comparison to constraints based on CH4 flux data alone. The improved model capability was attributed to the well-constrained parameters regulating the CH4 production and emission pathways. Our results suggest that the EBG modeling approach better characterizes CH4 emission and underlying mechanisms. Moreover, to achieve the best model results both CH4 flux and concentration data are required to constrain model parameterization